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Abstract 

Fractal and multi-fractal aspects of spatial intermittency in the energy cascade of electron 
magnetohydrodynamic (EMHD) turbulence is considered. Fractal and multi-fractal models 
for the energy dissipation field are used to determine intermittency corrections to the scal- 
ing behavior in the high-wavenumber (electron hydrodynamic limit) and low-wavenumber 
(magnetization limit) asymptotic regimes of the inertial range. Extrapolation of the multi- 
fractal scaling down to the dissipative microscales confirms in these asymptotic regimes a 
dissipative anomaly previously indicated by the numerical simulations of EMHD turbulence. 
Several basic features of the EMHD turbulent system which are universal and transcend the 
existence of the characteristic length scale d e (which is the electron skin depth) in the EMHD 
problem are highlighted. 
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1. Introduction 

The high-temperature plasmas in space (e.g. solar flares and magnetospheric substorms) 
and laboratory (tokamak discharges) have been found to be collisionless. An important 
aspect of a collisionless plasma is the enhancement by an order of magnitude of the magnetic 
reconnection rate (Yamada pQ). In situations where the spatial scales are shorter than 
the ion-inertial length di and time scales are shorter than the ion-cyclotron periodjl the 
ions do not have time to respond and merely provide a neutralizing background, and the 
dynamics are entirely controlled by electrons. A fluid description for the electrons then leads 
to the electron magnetohydrodynamic (EMHD) model (Kingsep et al. |3], Gordeev et al. 
[1]). The strongly sheared electron flows in the current sheets in EMHD undergo Kelvin- 
Helmholtz instability and lead to turbulence in EMHD (which is to be contrasted with 
turbulence generation/intensification via the tearing mode instability of current sheets in 
MHD). The energy cascade in EMHD turbulence proceeds directly even in two dimensions 
(2D), as in MHD turbulence, thanks to the Lorentz force on the electrons. Biskamp et 
al. [5], [6] did high resolution numerical simulation of decaying 2D isotropic homogeneous 
EMHD turbulence and found that the energy spectrum follows the Kolmogorov spectrum 
in the electron hydrodynamic limit (d e /£ ^> 1, d e being the electron inertial length) in spite 
of the fact that the whistler waves (which are generic to EMHD) would be expected to 
mediate the energy cascade. (A whistler-like relational implying an equipartition of energy 
between the poloidal and axial components of the magnetic field was however found to hold.) 
Celani et al. [8] further showed that a Kolmogorov 4/5th law type result also holds for the 
energy cascade in 2D EMHD turbulence. Numerical simulations of Boffetta [9] revealed 
the presence of spatial intermittency in EMHD turbulence - the energy dissipation field 
was found not to be uniformly distributed in space and the dissipative structures were of 
filament shape. Numerical simulations of Germaschewski and Grauer [TD] showed deviations 
from a Kolmogorov-type linear law of the characteristic scaling exponent of higher order 
structure functions further validating this aspect. Numerical simulations of Biskamp et al. 
[5] and [6] also showed that the energy dissipation rate in EMHD turbulence was apparently 
independent of the dissipation coefficients suggesting the possibility of a dissipative anomaly?] 
in the direct energy cascade in EMHD. 

In this paper, we consider fractal (Frisch et al. [11]) and multi- fractal (Frisch and Parisi 
[T2] ) models to describe the effects of spatial intermittency in 2D fully-developed EMHD 
turbulence. We will then extrapolate multi-fractal scaling in the inertial range down to 
the dissipative microscale and provide analytical evidence for a dissipative anomaly in the 
high-wavenumber (electron hydrodynamic limit) and low wavenumber (magnetization limit) 
asymptotic regimes. Several basic features of the EMHD turbulent system which are univer- 
sal and transcend the existence of the characteristic length scale d e in the EMHD problem 
are highlighted. 



2 In situ measurements in the solar wind have provided evidence of magnetic field fluctuations characterized by such spatial 
and time scales (Alexandrova et al. [2]). 

3 The extent of the whistlcrization in EMHD turbulence was numerically investigated by Dastgeer et al. [7|. 

4 The finiteness of the energy dissipation even in the limit the dissipation coefficients vanish constitutes a dissipative anomaly 
(persistence of symmetry breaking even in the limit the symmetry breaking factors vanish). There is experimental support 
(Sreenivasan I13p for this in 3D hydrodynamic turbulence. 
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2. Governing Equations of EMHD 



The 2D EMHD system of equations can be written in terms of two scalar potentials - the 
magnetic flux function A describing the magnetic field in the plane B = V x A i z and the 
stream function ip describing the electron flow velocity in the plane v e = V x ip i z , which is 
proportional to the in-plane current density (so ip also represents the out-of-plane magnetic 
field): 



the equation of generalized vorticity: 



4^ + !We-V)c--!-(B.V)J = ^V 2 u, (1) 



dt \ d\) ' m e n e c d^ 

• the generalized Ohm's law: 

1(A+ £j) + (v e ■ V) ^ + £jj = nV 2 A (2) 

where, 

-J = -V 2 A, u = -Vfy (3) 

c 

The number density n e is constant, in accordance with the incompressibility of the electron 
flow V • v e = which implies V ■ J = - this presupposes that the displacement current 
dE/dt is negligible. 

In the ideal limit {y and n =>- 0), equations (1) and (2) have the Hamiltonian integral 
invariant (upon appropriately non-dimensionalizing the various quantities (Biskamp et al. 
[5] and ©)) - 

H =\JJ K VA ) 2 + ^ + ^ + ( V ^) 2 }] dS W 

s 

S being the area occupied by the plasma. (4) shows that the magnetization effects introduce 
a characteristic length scale, namely d e in the EMHD problem. As a result, the latter exhibits 
some departures from the properties of MHD turbulence. One such feature is a decrease of 
the energy flux, leading to energy pileup of scales £ n ~ d e in the energy cascade. This could 
lead to an ordered quasi-crystalline phase signifying the appearance of long-range order in 
the system (similar to the case with geostrophic turbulence (Kukharin et al. [T4])). 

(4) implies, on noting a whistler-like relation!! ip ~ Aji holds between the poloidal and 
axial components of the magnetic field (Biskamp et al. [5] and [6]), that the energy per unit 
mass at length scale t is given by 

E~1?(l + %) (5) 



5 This relation also implies an equipartition in the energy contents of the in-plane magnetic field and velocity fluctuations 
(Alexandrova et al. [2]). 
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which, in the magnetization (d e jl <C 1) and the electron hydrodynamic (d e /£ 3> 1) asymp- 
totic regimes, leads to 

„ J d e /£ < i , R M 

(6a ' b) 

It is of interest to note that EMHD turbulence also exhibits some basic features which 
transcend the existence of the characteristic length scale d e in the EMHD problem. One 
such feature becomes apparent on applying the equilibrium statistical mechanics approach 
to the EMHD problem. 



3. Equilibrium Statistical Mechanics 

Consider an EMHD turbulence within a square which can be expanded into an infinite 
series of discrete wave vectors k n with stream function amplitudes ^/(k n ) related to each 
other via equations (1) and (2). The Fourier analysis of this system allows a formulation in 
terms of many degrees of freedom and hence leads to a consideration of this problem from 
the viewpoint of statistical mechanics. 

Application of equilibrium statistical mechanics to this system (Burgers [T5], Hopf |16j . 
Lee [T7] and Kraichnan [18]) requires the latter to be considered idealU This, in turn, 
requires a truncation in the Fourier space by retaining the Fourier modes lower than a cut- 
off wavenumber k max . This set of N wavenumbers conserves the energy (according to (5)), 
which is a quadratic rugged invariant, and survives the spectral truncation (because it is 
conserved by an interacting triad) 

\ E (1 + I 2 = canst. (7) 

If y n i (k n ) and y n2 (k n ) are the real and imaginary parts of each mode $ (k n ), the system 
can be represented by a point of m — 2N coordinates in a phase space and evolves ergodically 
in this phase space on the energy sphere, 



^ m 

-^2(l + k 2 a d 2 e ) y 2 a = const. 



Consider now a collection of such systems which is represented at each instant of time 
by a cluster of points in the phase space of density p(yi, ...,y m ,t). Since the total number 
of such systems and hence the volumes occupied by their representative points in the phase 
space are preserved, we have the Liouville Theorem: 

^,f^i-o ( 9 ) 



6 Formally, equilibrium statistical mechanics does not seem to be applicable to turbulence which, being dissipational, is in a 
non-equilibrium state. However, a turbulent system is believed to relax via nonlinear interactions toward equilibrium (which is 
confirmed by the numerical calculations of Orszag and Patterson |f 9|). Indeed, one may interpret the energy cascade to small 
length scales as a consequence of this tendency (Novikov |20| ~l . 
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Statistical mechanics seeks to explain the statistical behavior of a system in terms of 
its structural properties, such as the conservation of energy. This enables the equilibrium 
spectrum of EMHD turbulence to be predicted from the viewpoint of canonical ensemble 
averages. 

The elementary Gibbsian methods of statistical mechanics allow construction of equilib- 
rium solutions of Liouville's equation (9) as functions of the conserved quantities, such as 
the energy (8), via the Boltzmann-type distribution 

P (y u ..., y m ) = -e «=i (10) 

where a is a constant (interpretable as "inverse temperature" ) and Z is the partition function 
of the system, 

Z = ... e «=i dy 1 ....dy m . (11) 



The canonical ensemble average (p(yi, ...,y m ,t)) of an ensemble of the given system 

p (yi, y m , t) is stipulated to relax eventually toward this equilibrium distribution over the 

energy sphere (8) in the phase space. 

The mean variance of the mode a is given by 

(yi) = tf J ■■■ J vlp 2 ^ ^ + " ^ >Va dy 1 ...dy m (12a) 

or 

The energy spectrum is then given by 

E(k) ~ nk (1 + k 2 df) (|*(k)| 2 ) ~ 7ik, Vk. (13) 

(13) shows that EMHD turbulence, like the infrared regime of 2D hydrodynamic turbulence, 
exhibits the equipartition spectrumj^ E(k) ~ k, for small wavenumbers. This result signi- 
fies basic characteristic aspects of EMHD turbulence which transcend the existence of the 
characteristic length scale d e in the EMHD problem, as is also apparent in the following 
developments. 



4. Inert ial-Range Scaling Laws 

One may consider for the energy cascade in EMHD turbulence an inertial range of Kol- 
mogorov type which is in a state of statistical equilibrium and the energy is assumed to 
cascade smoothly through nonlinear processes in a stationary state. 

Consider a discrete sequence of scales, 

£ n ~£ -2- n ; n = 0,l,2,... . (14) 



7 The number of modes in 2D is proportional to 2irk. 
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Let us assume that we have a statistically stationary EMHD turbulence, where energy is 
introduced into the plasma at scales ~ £q, and is then transferred successively to scales 
~ £i, £2, • ■ ■ until some scale id is reached where dissipative effects are able to compete with 
nonlinear transfer. 

The energy per unit mass in the nth scale, according to (5), is given by 



f 2 

The rate of energy transfer per unit mass from the nth scale to the (n + l)th scale is 
given by 

e n ~ "7 W~ 1 + Tl) ( 16 ) 



'n 



where t n is a characteristic time of the nth scale, 

( " ~ T¥- (17) 

In the inertial range, we assume a stationary process in which the energy transfer rate is 
constant, 

e n = const = e, id < £ n < £o- (18) 



Using (18), (17) leads to 



e 



> 3 el> 3 / d?\" 1/3 



+ (19) 



Using (19), (16) gives 



from which we have 



//3 / ,2x1/3 



(20) 



„ J e 2 ^d e ^ in , d e /£ n <C 1 / Q1 „ 1 \ 

K ~\ eVHT,dJl n »l. (21a ' b) 

(21) leads to the following energy spectra (Biskamp et al. [5J and [5]), 

Efc ~ i e^k-^, kd e » 1. (22a ' b) 

The electron hydrodynamic limit corresponds to kd e ^> 1 while the magnetization limit 
corresponds to kd e <C 1. 



5. Spatial Intermittency 

The inertial range theory discussed in Section 4 does not take into account the spatial 
intermittency in EMHD turbulence that was revealed by the numerical simulations (Boffetta 
et al. [H] and Germaschewski and Grauer [ID])- Spatial intermittency effects would cause sys- 
tematic departures from the scaling laws (22) which use mean transfer rates. One may follow 
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Mandelbrot [21] and argue that the spatial intermittency effects in EMHD turbulence are 
related to the fractal aspects connected with the strongly convoluted dissipative structures 
(like the current sheets revealed in the numerical simulations [10]). This may be simulated 
in a first approximation by representing the dissipative structures via a homogeneous fractal 
with non-integer Hausdorff dimension D . This amounts to assuming the energy flux to be 
transferred to only a fixed fraction /3 of the eddies downstream in the cascade (Frisch et al. 

in])- 



5.1 Homogeneous Fractal Model 

Consider a discrete sequence of scales as in (15), but now assume that at the nth step, 
only a fraction /3 n of the total space has an appreciable excitation with a fractal dimension 

The energy per unit mass in the nth scale is given by 

/ d 2 



n 



En~P n K[l + f) (23) 



where, 

• (24) 

The energy transfer rate per unit mass from the nth scale to the (n + l)th scale is given 

by 

1 + | ■ (25) 

In the inertial range, the energy transfer rate is constant for a stationary process, so on 
using (24) and (25), and assuming the scaling behavior, 

v& n ~ C (26) 

we have, from (18), 



3a + 2- D 0(1) -2 = 0, d e /£ n < 1 
3a + 2- D 0(2) -4 = 0, d e /£ n > 1 

from which, 



(27a, b) 



3 ' 

Using (26) and (28), we have from (23), 



e 

(29) leads to the following energy spectra 



E(£) - { de J" , v ' de/in <<C (29a, b) 



E {k) ^^^ (30a,b) 

« e2 /3^5/3-l/3(2-A )(2) ) 5 fc4>>1 
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Observe that the intermittency corrections (2) < 2 J make the spectra steeper, as ex- 

pected. 

Noting that in the electron hydrodynamic limit (kd e ^> 1) the dissipative structures are 
typically vortex-filament like ^-Do {2) = OJ , and in the magnetization limit (kd e <C 1) they 

are typically current-sheet like \ Do m = lj ([10]), (30a, b) would lead to 



On the other hand, noting that the structure function of S p (£), order p, for the EMHD 
turbulence problem is defined in terms of the magnetic field in the magnetization limit 
(d e /£ <C 1) and the electron flow velocity in the electron hydrodynamic limit (d e /£ ^> 1), 

* p[q \(\S(d^/d£)(£)\P, d e /£^>l. [6 a > D) 

Using (26), and noting that the probability to belong to this fractal at scale £ goes like 
f- Do , (32) leads to 

°eo, d e /£^> 1. 

So, the characteristic exponent ( p is given by 



2 / 3 de 2/3 k~l kd e < 1 



S p (£) J_ 1)p+2 _ D > T,/L\ (33a, b) 



^-\{a-l)p + 2- D 0C2) , d e /£ » 1. ld4a ' Dj 



Using (28), (34) becomes 

2p 



(2-A, (1) ), 

C P = { ' (35a, b) 



P 
3 



(|-l) (2-D 0{W ), d e /£->\ 



which does not show a nonlinear dependence on p as is required of the characteristic exponent 
for large p. It is therefore necessary to consider the multi-fractal model [12J to address this 
issue. 



5.2 Multi- fractal Model 

Let us assume that the energy flux (or dissipation) is concentrated on a multi-fractal 
object Q12J) which is characterized by a continuous spectrum of scaling exponents a , a G 
/ = [a m i n , a max ). Each a & I has the support set S(a) C M 3 of fractal dimension /(a) such 
that, as £ =>- 0, the stream function increment has the scaling behavior - 

\8^(£)\~£ a . (36) 
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The sets S(a) are nested so that S (a') C S(a), for a' < a. The fractal dimension f(a) is 
obtained via a Legendre transformation of the scaling exponent of the pth order structure 
function of the electron flow velocity (or magnetic field), 

9 m ~ / / ^(«)^ p+2 " /(a) ~ , ^ « i r o 7a h N 

5pW t/d^a^^WM „ ^ (2 ), 4/£ » 1 ^ bj 

where the measure dfj,(a) gives the weight of different scaling exponents a, and £ 2 ~/( a ) 
represents the probability of encountering the set S(a) within a 2D circle of radius I. (37a, 
b) reflect the asymptotic scalings exhibited by (6a, b). 

One may use the method of steepest descent to extract the dominant terms in the integrals 
in (9), in the limit of very small £. This gives 

f a*p + 2 - /(«*), d e /£ < 1 

\(a*-l)p + 2 -/(«*), 4A>1 1 ' J 

where, 

/'(«*) = p. (38c) 

Next, in order to relate the singularity spectrum f(a) to the generalized fractal dimension 
(GFD) of the energy dissipation field, note that the energy transfer rate per unit mass at 
length scale £ is given by 

e( 4~7~\(d°/W^ (39a ' b) 

If the energy dissipation field is assumed to be a multi-fractal, the sums of the moments of 
the total energy dissipation U(£) ~ e(£)£ 2 occurring in N(£) squares of size £ covering the 
support of the measure e exhibit the following asymptotic scaling behavior (Halsey et al. 

[22]) 

V \U-UW ~ ~ I ^ ^(«)^" /(Q) , d e /£ « 1 f , 

2-, [Ul{ ^ 1 \jdfi(a)£^-^-^\ d e /£ » 1 (4Ua ' b) 

where D q is the GFD of the e-field (Hentschel and Proccacia [23]), arid we have assumed 
that the number of iso-a squares for which a takes on values between a and a + da is 
proportional to d(i(at)£~f( a ' . (40a, b) again reflect the asymptotic scalings exhibited by (6a, 
b). The dominant terms in the integrals in (40) may again be extracted, in the limit £ =>- 0, 
using the method of steepest descent, to give 

(n-l)D-S 3a*g-/(a*), d e /£ 1 (au b) 

where, 

/V) = 3g. (41c) 

The coincidence of the values of a* given by (38) and (41), for which the integrands in (37a, 
b) and (40a, b) become extremum, is insured by assuming a Kolmogorov refined similarity 
type hypothesis (Meneveau and Sreenivasan [21]) in the dissipative microscale regime. 
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Eliminating f(a) from (38) and (41), and putting q = p/3, we obtain 

{2-D p/3 ), 4/^«l 
^ = { ' (42a, b) 

| - (| - l) (2 - D p/3 ) , djt » 1. 

For a fractally homogeneous EMHD turbulence, 

(D d e /£<£l 
D ^ = \D^ V d e /i»l y P ( 43a > b ) 

(42a, b) reduce to 

(2-D 0(1) ), 4/^«l 
C P = <( * (44a, b) 



| - (f - l) (2--Do w ), 4/^»l 



3 V3 

in agreement with (35a, b). The energy per unit mass then shows the following scaling 
behavior, 

f e 2/3 rf -2/3 £ 4/3 + l/3(2-D ) rf /£ 
~ 2/3 4- 1/3^2 D ^ ' ^ 5a ' ^ 

and the energy spectra are, 

f e V3 4 -/3^/3-l/3(2-A, (1) ) 5 fc4<<1 
I ^3^/3-1/3(2-^^ fc4>>1 

in agreement with (30a, b). 

5.3 Multi-fractal Scaling at the Dissipative Microscale 

We now consider extrapolation of the multi-fractal scaling in the inertial range discussed 
in Section 5.2 down to the dissipative microscale by assuming that an inertial behavior 
persists at scales smaller than d e - this assumption may be justifiable for tenuous plasmas 
like those in space (d e ~ 10 km for the magnetospheric plasma). 

On taking into account the spatial intermittent character of the energy dissipation field, 
the dissipative microscales £d (1) {2) , 

n^d" 1 

Cd w ~ £l/2 e , d e /£ < 1 (47a) 

3/4 

^W~^74. de/i^l (47b) 

(along the lines of the development of Paladin and Vulpiani [25] and Nelkin [26J for the 
hydrodynamic case), are found to exhibit the scaling behavior, 

~ R~ 1/a , d e /£ « 1 (48a) 
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Co m ~ R^ /a , d e /i » 1 (48b) 
where R m and Rh are, respectively, mean magnetic and hydrodynamic Reynolds numbers, 

^m^.^w&l (49 ) 

rj v 

e is the mean energy dissipation rate, and rj is the resistivity and v is the kinematic viscosity 
of the plasma. (Observe that the mean magnetic and hydrodynamic Reynolds numbers are 
both dependent on the electron skin depth d e .) The identity of the scaling exponents in 
the two opposite asymptotic regimes is symptomatic of certain universal features in these 
regimes, as seen further in the following. 

The moments of the electron-flow velocity (or magnetic field)-gradient distribution, 

Ap - 1 (\d 2 iP/dx 2 \P), d e /£ > 1 lbUa ' bj 



are then given by 



, (/^ m )- i[( "- 1, '' +2 - /( ' >,, , dje«i 

' \fdKa)(R h r i[(a - 2, » +2 - f,a,] ,d c /e>>i. 



Ap 



In the limit of large R m and Rh, the dominant exponents in (51) correspond to 

a* [p - f(a*)] = (a* -l)p + 2- f(a*), d e /£ < 1 (52a) 

a* [p - f(a*)] = (a* -2)p + 2- /(a*), d e /l > 1. (52b) 

The coincidence of the values of a* given by (41) and (52) for which the integrands in 
(40a, b) and (51a, b) become extremum, is again insured by assuming the Kolmogorov 
refined similarity type hypothesis ([21]) in the dissipative microscale regime. (52a, b), in 
conjunction with (41a, b), lead to 

D Q (p-3)-3p + 6 

(Rm) Dq , where Q = ° Q + J " 2 , d e /£ « 1 

D Q (p-3) -6p + 6 Q (53a ' b) 

{Rh) D Q + 2 , where Q = Dq ^ 2P ~ 2 , d e /i » 1 

Dq + I 

So, the mean energy dissipation has the following scaling behavior, 

V A 2 ~ (Rm)°, d e /£ < 1 (55a) 

vA 2 ~ (R h )°, d e /£^>l. (55b) 

(55a, b) implies an inviscid dissipation of energy in the electron hydrodynamic limit and a 
non-resistive dissipation of energy in the magnetization limit and hence a dissipative anomaly 



from which, 
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in high- and low-wavenumber asymptotic regimes of EMHD turbulence in confirmity with 
DNS ([5], [6]). Note further from (53a, b) that the energy dissipation field in these asymptotic 
regimes has the GFD Dq equal to the information entropy dimension D\. The dissipative 
anomaly signifies another basic characteristic aspect of EMHD turbulence which transcends 
the existence of the characteristic length d e in the EMHD problem. 



6. Discussion 

One may view the energy dissipation rate e to be the order parameter a la Landau [27] for 
the EMHD turbulence problem because it appears to indicate the degree of broken symmetry 
and exhibits fluctuations in the presence of spatial intermittencyJl Further, noting that the 
critical point for EMHD turbulence corresponds to the limit R m and Rh =>■ oo, the non-zero 
limiting value of e, as the critical point is approached appears to validate this view. Indeed, 
one may define the critical exponent o|l (Shivamoggi [28]) for this problem by 



(R m ) a , R m oo, djl < 1 
(R h ) , Rh => oo, d e /£^> 1 



(56) 



where, as per (53a, b), 

a = 3(Q - 1), W e /£. (57) 

Comparison of the above results (see Table 1) with the corresponding results for various 
FDT systems in fluid and plasma dynamics (Shivamoggi (33] - [35]) indicates that the energy 
(or enstrophy in 2D hydrodynamic FDT) dissipation rate e is the right choice for the order 
parameter for the FDT problem^ with an apparently universal form for the critical exponent 



8 This perspective therefore allows (Shivamoggi 1281 ) the Kolmogorov type inertial range theory described in Section 4, which 
assumes e to be uniform, to be appropriately regarded as a mean field theory d la Landau |27| . 

9 One of the goals of critical phenomena formulation of the turbulence problem (Nelkin I29| . Yakhot and Orszag 30 , Eyink 
and Goldenfeld 1311 . Esser and Grossmann 32 , Shivamoggi 28 ) has been to determine the critical exponents that are intrinsic 
features of the turbulence dynamics and are not artifacts of the large-scale stirring mechanisms. 

10 It may be mentioned that different choices (Nelkin [29) . Rose and Sulem [36]) have been considered for the order parameter 
for the FDT problem; the present choice seems to be appealing because it agrees with all the implications posited in Landau's 
order parameter concept (see also footnote 8, as well as remark above equation (56)). 
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a given b}q 

o- = a (Q - 1). (58) 

The variations in the amplitude a reflect the residual effect of variant cascade physics in the 
diverse FDT systems^! Observe in Table 1 that the energy (or enstrophy) dissipation fields 
in the various FDT systems have the GFD Dq equal to the information entropy dimension 
D\ (because, corresponding to p = 2, the GFD index Q turns out to be unity for all these 
FDT cases). 

Further insight can be gained into this aspect by looking at the probability distribution 
function (PDF) of the electron- flow velocity (or magnetic field) gradient. In order to derive 
the PDF of the electron-flow velocity (or magnetic field) gradient, note that the scaling 
behavior of the dissipative microscales, on using (47a, b), is given by 

^-j , d e /£ « 1 (59a) 

j\ , deft » 1 (59b) 

where ipo is the stream function increment on a macroscopic length L. 

The scaling behavior of the electron-flow velocity (or magnetic field) gradient is then 

JL~^ a ,(«-D/«, dji«i 

(6 ° a ' b) 

^(2) 

The PDF of the electron-flow velocity (or magnetic field) gradient may then be determined 
in terms of that for the characteristic stream function increment ipo for large scales as follows, 

P(a;a) = P(V>b)^. (61) 
as 



11 The critical exponent cr may be connected with the critical exponents 7, v and rj introduced by Rose and Sulem 1361 
according to 



as follows, 
with 

Noting, from (28) and (29), that 
and 

we obtain 



S 2 (£) - tn ~ i<* so r, = C P 
a = 3(7-1) 
Q= 7 = f(2-J7). 

_ J3/2 , d e /e < 1 
v \3/4 , d e /e > 1 

_ ("4/3 , d e /e < 1 

v ~ \2/3 , de/i^l 



Q = 7 = l 
a = 



v d e /e 



as required. 

12 This is totally in accord with the idea of universality which implies that near a critical point all systems can be grouped 
into a relatively small number of classes (depending on the specific dynamics) with identical critical exponents within each class 
(Hohenberg and Halperin |37] ). 
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FDT case 


Critical Exponent a 


Generalized Fractal Dimension Index 


3D incompressible FDT 


3(0-1) 


D Q + 2 P -3 
W Dq + 1 


2D incompressible 
FDT-enstrophy cascade 


3(0-1) 


D Q + 3 P -2 
* Dq+A 


3D compressible FDT 




Q 


27 

Dq + —^rp - 3 
7+1 

47 

Dq + - 3 
7 + 1 


3D MHD FDT 


2(0-1) 


2D Q + 3p-6 


2D EMHD FDT 


3(0-1) 


0= < 





Table 1: Critical exponents for various FDT cases. 
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Taking Pijpo) to be Gaussian, 
and using (60a, b), (61) leads to 



~ e 



-Vg/2<Vg> 



P(s: a) ~ < 



l-Q 



V 



2(l-a)| 1 2a 



11/2 



2-a 



2 < ^ > _ 

■ u 2(2-a)\ s \a 



d e /£ < 1 



4/^ > l. 



For EMHD turbulence, on noting from (26) and (28) 



a 



2/3, d e l I < 1 

4/3, c? e /£ > 1 



(63a, b) become 



Pis) ~ < 



1/3 



2/3 I „ I 4/3 



2\ V3 



2 < % > 
^4/3^14/3 ' 



2 < r > 



d e l I < 1 



4/^ > 1. 



(62) 



(63a, b) 



(64a, b) 



(65a, b) 



The identity of the \s (-dependence exhibited by P(s), as per (65a, b), (which is also the 
same as the PDF for the velocity gradient for 3D hydrodynamic turbulence given by Frisch 
and She [38]), appears to be consistent with the demonstration of dissipative anomaly, as 
per (55a, b), in the asymptotic regimes (this is validated further by the critical exponent 
(57) for EMHD). The stretched exponential decay of the PDF exhibited by (65) has also 
been indicated by the numerical simulations ([9] and [TO]). 
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